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Abstract — Finite Element Analysis (FEA) using Finite Element Method (FEM) was 
developed over 70 years to solve the complex elasticity and structural analysis problem in 
civil and aeronautical engineering. Application of FEA is being expanded to simulation in 
electrical engineering also to solve the complex design problems. The circuit theory models 
for designing transformers are not much accurate in determining the transformer 
parameters such as winding impedance, leakage inductance, hot spot temperature etc. The 
physical realization of these parameters is needed on a prototype unit. The finite element 
method can play a vital role in deriving these parameters without any physical verification. 
An effort has been made in this paper to show the effectiveness of finite element method in 
determining the above said parameters while designing the transformers - both oil cooled as 
well as dry type - for power and distribution sectors as well as to analyze and detect the 
internal faults in the transformer. 

Index Terms — Finite Element Method, Finite Element Analysis, Quasi-static Finite Element, 
Leakage Inductance, 

L Introduction 

The finite element method (FEM) rapidly grew as the most useful numerical analysis tool for engineers and 
applied mathematicians because of its natural benefits over prior approaches. The main advantages of FEM 
are, it can be applied to arbitrary shapes in any number of dimensions, the material properties can be non- 
homogeneous (depend on location) and/or anisotropic (depend on direction). The way that the shape is 
supported (also called fixtures or restraints) can be quite general, as can be the applied sources (forces, 
pressures, heat, flux, etc.). The FEM provides a standard process for converting the governing energy 
principles or governing differential equations in to a system of matrix equations to be solved for an 
approximate solution. For linear problems such solutions can be very accurate and quickly obtainable. 
The objective of applying finite element method while designing a transformer is, to estimate the temperature 
rise and hot spots under non linear load conditions, to determine the impedance value and short circuit forces 
to identify the electrical stresses and weakest insulation points and to estimate the correct insulation level to 
avoid the HV failures. 

This paper discusses the use and effectiveness of finite element method in designing the power transformer as 
compared to magnetic circuit theory [12], internal winding fault detection and analysis [3, 4], a three 
dimensional finite element analysis of electric fields at winding ends of dry type transformer [9], hot spot and 
life evaluation of power transformer [6, 13], leakage inductance calculations [7] and transformer over heating 
under non linear load conditions [8] 
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II. Transformer Design using fem 



In designing a transformer, especially power transformers, parameters like impedance of the transformer, 
temperature rise and hot spot evaluation, rated short time overloading capacity and adequate insulation level 
are of prime importance. Determination of these parameters at design stage has a very big impact on the cost 
effectiveness of the transformer. 

In a paper presented by S C Bell and S Bodger [2, 15], a comparison has been made between magnetic circuit 
theory and finite element method for designing the power transformer. It summarizes the reverse method of 
transformer design. 

In this paper, first the models for the resistive and inductive reactance components of the Steinmetz 'exact' 
transformer equivalent circuit have been developed from fundamental theory as presented in [13]. Then two 
and three dimensional linear and non linear magneto static finite element models were introduced as an 
alternative model for the inductive reactance components. 

To demonstrate the reverse design method, the author have designed, built and tested two single phase, 50 
Hz, high voltage transformers, the results of which are summarized in table I. 

Table I. Calculated And Measured Equivalent Circuit Parameters For Sample Transformers 



Method 


Equivalent Circuit Parameters in Q 




Rc ! ! Xm \ %l 


Transformer 1 

Primary 240 V, Secondary 6240 V, 200 VA 


DM 


3388 


10.0 


1987 


2.8 


CTM 


1342 


11.5 


1383 


1.9 


LFEM 






1905 


1.6 


NLFEM 






1883 




Transformer 2 

Primary 14 V, Secondary 4560 V, 617 VA 


DM 


18 


0.043 


41 


0.012 


CTM 


9.9 


0.055 


20 


0.016 


LFEM 






25 


0.015 


NLFEM 
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R c = Transformer core resistance 

R w = Transformer winding resistance 

X m = Magnetizing reactance of the transformer 

X l = Leakage reactance of the transformer 

DM = Direct Measurement 

CTM = Circuit theory model 

LFEM = Linear finite element model 

NLFEM = Non-linear finite element model 

The results show that the non linear finite element model most accurately calculated the magnetizing 
reactance value of the two sample transformers. For transformer 1 the non linear model seems to be less 
accurate however this may be due to the approximations made in the geometry of the finite element model. 
D. Azizian, M. Vakilian, and J. Faiz [7] have introduced analytical and FEM based models for 
electromagnetic modeling and inductance calculation in dry type multi-winding traction transformers. The 
accuracy of these models is verified against the experimental data for a traction transformer having 
specifications as given in table H 

The test results are as shown in table III from which it is evident that the FEM has better and most accurate 
results. The accuracy of axi-symmetric 2D FEM model goes on increasing as we increase the complexity of 
modeling from simplest core (2D-S) to full core (2D-F) to improved core (2D-I). The accuracy of 3D model 
is much higher than the 2D models. 

The experimental results of this paper also yield that the proper modeling of the transformer core further 
helps in determining the value of leakage inductance due to change in the cross section area of the core. This 
is especially very helpful in designing the converter duty transformers where leakage inductance between the 
transformer windings is of utmost importance to decide the capacitance value for the LCL filter. The leakage 
inductance of the transformer itself forms a part of the filter circuit. 
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Table II. Specifications Of Typical Traction Transformer 



Power Rating 


Vector Group 


4000 KVA 


DdO/yll 


HV Line Voltage 


HV Line Current 


(H1||H2) = 20 KV 


(H1||H2)= 115.5 A 


LV Line Voltage 


LV Line Current 


(LI || L2) = 750 V 


(LI || L2)= 1540 A 



Table III. Leakage Inductances Of Traction Transformer 



Test/ 


Leakage Inductances (%) 


Measurement 


Ll- 


Ll- 


Ll- 


L2- 


L2- 


Hl- 


Method 


Hl 


H2 


L2 


Hl 


H2 


H2 


Analytical 
Results 


5.9 


72.5 


55.2 


71.7 


88.1 


6.0 


2D-Simplest 
core 


5.679 


44 


45 


44 


42 


5.686 


2D-Full Core 


5.808 


111 


112 


112 


110 


5.84 


2D- Improved 
Core 


5.794 


98 


100 


99 


97 


5.834 


3D Model 


5.892 


88 


88 


88 


88 


5.914 


Experimental 


5.797 


74 


73 


72 


72 


5.943 



III. Life evaluation of Transformer using fem 



One of the most important parameters governing the life expectancy of a transformer is the hot-spot 
temperature value. Stray losses and non linear loads in a transformer are one of the main contributing factors 
in creating such hot spots which all together decide the life of a transformer. 

The stray losses in a transformer are caused by the time variable leakage flux which induces emf and 
circulates eddy currents in the winding as well as conducting parts of the transformers such as clamps, core, 
tank wall etc. Evaluation of stray losses can be done quite accurately by FEM as discussed by A. S. Reddy 
and M. Vijaykumar [6]. 

Dejan Susa [13] has developed the models to determine the hottest spot in a transformer based on heat 
transfer theory, application of the lumped capacitance method, the thermal-electrical analogy and a new 
definition of nonlinear thermal resistances at different locations within a power transformer. The changes in 
oil viscosity, loss variation with temperature and changes in transformer time constants due to changes in oil 
viscosity were also accounted for in the thermal models. The results showed that the top oil temperature time 
constant is shorter than the time constant suggested by the present IEC loading guide, especially in cases 
where the oil is guided through the windings in a zigzag pattern for the ONAN and ONAF cooling modes. 
The models are validated using experimental results, which have been obtained from a series of thermal tests 
performed on a range of power transformers. 

With development of the electronic equipments such as Computers, UPS, and High Frequency Drives for 
motor loads, arc furnaces, Electronic Ballasts, Compact Fluorescent Lamps, etc. the harmonic content in the 
power distribution network has increased tremendously. These highly non linear loads are present not only in 
the industrial sectors but also in the commercial sector and the first victim to any such load is always a 
transformer feeding them. These non linear loads result into the overheating of the transformers. The non 
linearity of the loads is best judged by the K factor [1] as described by equation (1) and accordingly the 
transformer required will be designated as K-rated transformer. 

The K-rated transformer does not mitigate the harmonic but is capable enough to sustain the overheating due 
to such non linear loads. The design and application considerations along with testing approach based on UL, 
NEMA and IEEE standard C57. 110-1986 are best described by L. W. Pierce, member IEEE [5]. 

K = £(yfV (1) 

~ Vrms/ 

Where, 

h = order of harmonic 
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l h = RMS value of the current for h order harmonic 
I rms = RMS value of the total load current 

A three dimensional finite element method using a magnetic scalar potential formulation is applied [8] to 
compute the magnetic field in free and iron spaces. The calculation is then combined with a mixed analytical 
and numerical form of the electrical circuit equation to take into account the skin and proximity effects in the 
rectangular windings in dry type transformers under non linear load conditions. A two step FEM using the 
both reduced and total magnetic scalar potentials throughout a penalty term [10] is proposed in [8]. 
Accordingly the magnetic field into the air gap (Qa) and in the core (Qc) based on the reduced and total 
magnetic scalar potential <|) follows the governing equations: 

lnQ. a :H = Hj — grad(0)andlnQ c : H = —grad(0) 

(2) 

Where, the source magnetic field Hjcan be computed using either the Biot-Savart law or a fictitious field for 
the known current distribution [11]. The relationship given in equation (2) has been solved from a two step 
formulation using a difference magnetic field h= H - H , considered as being a disturbance of H due to the 
core saturation [10]. In first step the magnetic permeability \x is assumed to be infinite and a reduced 
magnetic scalar potential is used to calculate the magnetic field H . In the second step, a total magnetic 
scalar potential formulation is applied with Ho as the new magnetic source field in order to calculate the 
difference field h. 

Based on the magnetic field described above and the empirical formula of the ac resistance component as in 
[14], the additional losses in the winding due to the harmonic currents have been modeled in [8] on a sample 
10 KVA dry type distribution transformer. The nominal current of the transformer is applied to the windings 
for all frequencies in the above paper. However with non linear loads the harmonics do not own the same 
amplitude. The power losses are then need to be weighted and modeled in accordance with the current to 
estimate the exact value of temperature rise. 

The other major factor that decides the life of a transformer is the designing of insulation system. The 
determination of insulation level can best be identified by determining the electric fields at the winding ends. 
Using three dimensional FEM and ANSYS software, J. Hong, L Heyun and X. Zihong [9] have estimated the 
electric fields at the winding end in 10KV SG10 dry type transformer. The similar approach can be extended 
to other higher rating transformer too at the design stage itself to estimate the correct insulation level and 
consequently increase the life of a transformer. 

IV. Fault detection in Transformer using fem 

Majority times the faults in a transformer are turn to turn, disk to disk or turn/disk to earth short circuit. This 
occurs mainly because of the aging of the insulation or displacement of the insulation during transportation or 
maintenance. This leads to the overheating of the transformer and finally results into the transformer failure. 
Study shows that around 70-80% of the transformer failures are due to the short circuit between turns. 
Considering this damage occurring in the transformer and time and cost involved in rectifying the same, it 
seems that simulation involving modeling of the transformer for detecting the fault is the most economical 
and convenient way. 

One of such methods of modeling distribution transformer with internal short circuit faults using Finite 
Element Analysis (FEA) is presented in [4]. Based on the physical information of the transformer, the finite 
element model for a normal transformer or a transformer with an internal fault was implemented by 
commercially available software. The resulting circuit model was exported and used in a circuit analysis 
package to study the terminal behavior of the transformer. The observations of the simulation are compared 
with the experimental results and it shows that the FEA model can provide an accurate estimation of the 
internal winding faults. 

The electromagnetic quasi-static finite element approach [3] describes the detection of winding short circuit 
faults with the help of frequency response analysis. The principle of coupled circuit approach and 
electromagnetic quasi-static analysis approach were applied along with FEA for determining the short circuit 
current under different fault conditions. 

A case study on 30 MVA, 63 KV/20 KV, Yndl, and 50 Hz transformer has been discussed in this paper [3]. 
Parameters of transformer are estimated by means of finite element analysis and utilized in circuit based 
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model and the input impedance is calculated in wide band frequency. In addition to classifying and analyzing 
main types of short circuit according to IEEE standard C57. 140, frequency response of disk-disk short circuit 
state is also investigated in several points in the HV winding to identify the location of the short circuit fault 
along the winding. The results have been summarized in terms of deviation of resonance component which 
showed that the first resonance of the input impedance due to short circuit moves orderly giving a better idea 
about the location of the short in a winding. 

V. Conclusion 

The discussion above reveals that the finite element method is an efficient tool for designing the transformer. 
Whereas the determination of impedance, winding hot spots temperature, insulation class and insulation level 
etc. are the major challenges while designing a transformer, the FEA can serve a vital role and may provide a 
cost effective and efficient solution to the problem. The preliminary results obtained from finite element 
model may be useful to develop cost effective and efficient design. However from all above discussions, it is 
evident that a three dimensional finite element model is predominant over a two dimensional model. While 
initiating the analysis a two dimensional model may be adopted but for the betterment of results the analysis 
must be completed by developing a three dimensional model only. 
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